A Three Dimensional Lattice of Ion Traps 
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We propose an ion trap configuration such that individual traps can be stacked together in a 
three dimensional simple cubic arrangement. The isolated trap as well as the extended array of ion 
traps are characterized for different locations in the lattice, illustrating the robustness of the lattice 
of traps concept. Ease in the addressing of ions at each lattice site, individually or simultaneously, 
makes this system naturally suitable for a number of experiments. Application of this trap to 
precision spectroscopy, quantum information processing and the study of few particle interacting 
system are discussed. 
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The trapping of ions with static [l| and time varying [2| 
fields has been central to many important developments 
in physics and chemistry for the past 50 years. Ion traps 
have been applied to precision spectroscopy ll| , preci- 
sion mass determination mass analysis U atomic 
ion clocks |5| and most recently in quantum information 
processing [6| . Ion trap experiments most typically focus 
on the ability to trap (1) a single ultra-cold ion for ex- 
tended interrogation, (2) a collection of ions in the form 
of a spatially extended cloud 0] or (3) an ordered crystal 
of ions d, 0] . Multiple species of ions can also be trapped 
simultaneously [l(| • Common to all of these experiments 
is the study of ion(s) within isolated ion traps. Recently 
experiments have been performed to communicate be- 
tween trapped ions in remote experiments 11| . Different 
trap configurations have been developed to meet specific 
requirements, starting from conventional 3-dimcnsional 
(3-D) traps with hyperbolic electrodes Q to cylindrical 
traps fl2j l. linear traps [l3| and various kinds of planar 
traps [li, [H, HE 03 ■ 



In this Letter we propose a novel configuration for 
ion trapping which is three dimensionally symmetric by 
its geometry. It resembles the configuration used by 
Wuerker et al. [llj in 1959. An important consequence 
of the present structure is its suitability for the construc- 
tion of several compact, stacked, virtually independent 
ion traps in the the same apparatus. Here we discuss the 
fundamental structure underlying the proposed trap and 
its operating scheme. Parameters for both the isolated 
fundamental trap and the system of stacked traps are 
determined and compared, establishing the robustness of 
the configuration. The motion of ions in the traps is char- 
acterized, building a strong case for several experimental 
applications with this trapping scheme, some of which 
are sketched to illustrate its versatility. 

The fundamental building block for the Lattice of Ion 
Traps (LIT) is illustrated in Fig. [Tfa). The trap is con- 
structed with 3 parallel pairs of cylindrical electrodes 
along the cartesian axes. When extended this creates 
compact, well isolated individual traps with large optical 



access, permitting optical communication between differ- 
ent traps and ease of construction. Typical dimensions 
could be electrode diameter 0.5mm and a center-to-center 
separation of 5mm. These are the dimensions used in 
our simulations. The dimensions can be changed with- 
out qualitatively changing the characterization below. 

Time varying potentials create the 3-D ion trap, 
when applied to the equipotential, parallel electrodes in 
Fig.QJa). The phase relation for the amplitude of the ap- 
plied radio frequency (RF) field between electrode sets 
along x, y and z is illustrated in Fig. [TJb). This field 
configuration creates the rotating saddle potential which 
traps the ions in a 3-D symmetric configuration over one 
full cycle of the RF. This symmetry lends itself to exten- 
sion, to form the LIT illustrated in Fig. [T{c). In general 
any I x m x n number of equivalent traps can be con- 
structed by building an appropriate number of electrodes. 
The RF phase relation remains identical to the funda- 
mental trap for the extended sets of parallel electrodes. 
Fig. QJc) illustrates the 3x3x3 version of the LIT as it 
is the lowest order cubic arrangement that provides lat- 
tice sites with different near neighbor coordination. The 
four distinct locations of ion traps in this arrangement 
are Body Center (BC), Face Center (FC), Edge Center 
(EC) and Corner (C). We show below that the site spe- 
cific effects on the trapped ions in individual traps are 
too small to be significant. 

The ion trap operation is simulated for the 40 Ca + ion. 
The potentials for the ion trap(s) are generated using 
SIMION. The equations of motion in the time varying po- 
tentials are solved numerically using Mathematica. The 
rotating saddle potential results in ion motion with a 
macromotion frequency v m whose amplitude is modu- 
lated by the driving frequency of the trap, v r f. Since v m 
can be significantly less than v r f care was taken in solving 
for the trajectory of the ion in three dimensions. Char- 
acterization simulations for the isolated ion trap and the 
LIT potential were done by solving for the full 3x3x3 
cell configuration. The isolated trap solutions are then 
compared with the corresponding solutions within the 
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FIG. 1: (color online) Structure and operation of the 3-D ion 
trap. Part (a) shows the basic electrode configuration of the 
ion trap. Three pairs of parallel electrodes are arranged along 
the cartesian axes and are electrically contacted. The relative 
phase difference of 27r/3 between the orthogonal electrode sets 
is shown in (b), where (Vi,V B ,V z ) are represented by (thin, 
thick, dashed) lines respectively. The extension of the trap in 
(a) to 27 traps in simple cubic arrangement is illustrated in 
(c). 4 distinct types of coordination symmetries for traps at 
different locations emerge (d). Representative cells of different 
coordination are identified as BC(6), FC(5), EC (4) and C(3), 
where the bracketed number is the location dependent nearest 
neighbor cells. 
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FIG. 2: (color online) Details of trap characterization. Part 
(a) illustrates the region of stability in the |vo| vs. V r / 
plane at v r j — 1MHz. The region contained within the the 
curve is the stable region. The dashed(black), thin(red) and 
thick(grey) lines are the stability bounds for the single trap of 
Fig- da), the most configurationally distinct BC and C cells 
in the LIT respectively. Part (b) plots the maximum initial 
velocity that is trapped | vomeral and (c) the maximum ampli- 
tude of the RF field which traps the 40 Ca + ion as a function 
of v r f. The line is the least square fit to the points, for (b) 
av r f and j3v^f for (c). 



individual traps at the BC, FC, EC and C locations. 

Operationally, apart from the dimensions of the trap, 
3 parameters characterize the symmetric 3-D ion trap. 
These are frequency of operation, amplitude of the RF 
and the trap depth for the ion. Equal amplitudes of 
the RF voltage V, / applied as in Fig. QJb) yields v m ^ x = 
v m ;y = ^m ; z- No dc field is applied to the electrode in our 
trap configuration. As a consequence of this the standard 
Paul trap representation in terms of stability parame- 
ters a and q, related to the Mathieu differential equation, 
which governs the motion of ions in a RF trap [2| is not 
relevant as a — 0. However 

q = 4QV rf /M{2irv rf ) 2 R 2 Q , 

where Q represents the charge state of the ion, M the 
mass and Rq = 0.0025m the distance of the electrode 
from the trap center i.e. the size of the trap, remains a 
good expression for trap characterization. 

The nature of the present trap is characterized by re- 
lating the RF amplitude V r /, the magnitude of the ini- 
tial velocity of the trapped ion |v | and v r f. Fig. 
illustrates the boundary of the region of stability in the 
vo| — V r / plane, where u r f = 1MHz, for the single trap, 
BC and C locations. The initial position of the ion is at 



the trap center and initial velocity is Vrj with equal com- 
ponents in the 3 orthogonal directions. The small varia- 
tions in the trapping region in Fig.[2]Ja) arise due to trap- 
ping potential at the edge of the physical trap. A barely 
trapped ion samples the outer reaches of the individual 
trap potential where, in addition to intrinsic anharmonic- 
ity, the trap at its edges, is susceptible to its immediate 
neighborhood. Thus the differences due to the trap lo- 
cation in the LIT are minor. The top axis of Fig. [Ha) 
marks out the macromotion frequency fraction 
of v r f , for the isolated trap and is only dependent on the 
magnitude of V r / . As V r / increases the amplitude of the 
micromotion increases rapidly, thus taking the ions to the 
trap boundary and ejecting them. Deeper trap potentials 
therefore do not imply higher initial kinetic energies of 
trappable ions, as illustrated by the fall in the |vo ;ma x| 
for higher V r / values. Thus the suitable operating point 
of the ion trap is the initial region of linear response of 
|vo;ma:c| t° V r / with q < 0.29 for V,-/ < 8V at v r f = 
1MHz. Generating multiple |vo| — V r f stability plots for 
different v r f and compiling the results provides the data 
for Fig. [2]Jb) and (c). Fig.[2[b) illustrates that |v ; maa :| 
scales as v r f , quantifying the maximum trappable initial 
velocity of the ion. Fig. \T[.c) demonstrates that the max- 
imum trapping voltage V r y ;ma;r scales as v 2 ? consistent 
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FIG. 3: Power spectrum of the x component of the trajectory 
for a low | vo| ion computed at v T f — 1MHz, V r f = 3V. 
A strong macromotion peak at v m and significantly weaker 
peaks at u r f, v r f ± I'm, 2v r f and 2i/ r f ± v m are seen. 
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FIG. 4: (color online). Power spectrum of the x component 
of the trajectory for a low |vo| ion computed at v r j — 1MHz, 
V r f = 3V for the ion trap with split electrodes (inset). The 
resulting spectrum is qualitatively similar to that of Fig. [3] 



with the behavior for constant q, here q — 0.908, which 
barely traps the zero velocity ion. Thus the bounds for 
the operation of the present ion trap are well character- 
ized via the scaling functions shown. 

Symmetric operation of the trap gives identical macro- 
motion frequency v m , in the orthogonal x, y and z direc- 
tion. The micromotion of the trajectory components is 
synchronous with the drive field frequency v r j. At the 
trap center and for a large range of |vn|, v m is unchanged 
and the power spectrum Fig. [3] remains free of sidebands, 
indicating that the trap is harmonic. Due to the large 
difference between v m and v r f in the q < 0.29 range of 
operation, it is reasonable to separate the ion motion into 
2 components and the slow motion at v m can be consid- 
ered separate, while averaging over the fast oscillation at 
v r f. This is indeed exhibited clearly in the amplitude of 
the power spectrum of the fourier transform along the 
representative x direction in Fig. [3J In the small Vo and 
low q region of operation, v m is the same across traps 
with different coordination numbers. 

The error in v m determination from our simulations is 
< ±0.01i<v/- The likely source of this error is the nu- 
merical solution to the Laplace equation, which at our 
computational convergence limit is marginally different 
from the true value of the potential. The complexity of 
the electrode structure here does not allow a fruitful com- 
parison for the motion of a trapped ion in the computed 
trap and its idealised analytic form. 

Structurally the electrodes of the trap illustrated in 
Fig. [TJa) and (c) do not allow a face on view of the ion. 
This can be easily circumvented by splitting each of the 
six electrodes into two closely spaced electrodes as shown 
in Fig. IDJa). The resulting trap is similar in performance 
to the trap in Fig. QTa) as illustrated in Fig. |4jb) by the 
power spectrum of the x component of the ion trajectory. 

Loading the LIT with the desired number of ions is 
most efficiently achieved by resonance enhanced two pho- 
ton ionization. Briefly, the RF fields could be applied 
while a vapor of the parent neutral floods the lattice. Ap- 
plication of the two laser frequencies participating in the 



ionization, along two distinct directions of easy viewing 
through the lattice would create ions only at the beam in- 
tersection. Thus ion(s) can be loaded in a controlled way 
in each ion trap cell, close to the center of the individual 
traps. Laser cooling of ions to ultracold temperatures, 
now well established, is easily done along the easy view- 
ing directions. Below we explore some applications which 
exploit the symmetry and multiplicity of the ion trap. 

The most precise spectroscopic measurement on a sys- 
tem would be on an isolated, trapped single particle, 
nearly at rest. However in the specific case of a single 
ion confined within an ion trap, a serious constraint is 
presented by the low flux of the photons and hence the 
signal/noise (S/N) in the detection. Adding more than 
one ion per ion trap leaves the system vulnerable to per- 
turbations that could compromise the precision of the 
measured transition frequencies. With the experimen- 
tal setup proposed here, we have the ability to preserve 
the isolation of the individual ions in each trap, while 
increasing the detected fluorescence rate. 

Specifically the performance of single ion optical clocks 
may be significantly improved when a single ultra-cold 
ion is trapped within each cell of the LIT. The stability 
of clocks expressed by the Allen variance 

0= (Aiy/iy)(S/N)-^ 2 x t" 1 / 2 , 

where v is the transition frequency and Av the observed 
linewidth and r is the integration time. In the arrange- 
ment illustrated in Fig.QJc), 27 such trapped ions can be 
interrogated simultaneously, increasing the S/N. This 
results in a decrease in the short term instability by about 
a factor of 5. 

Information processing with ion traps allows the most 
precise control over qubits thus far. However one practi- 
cal limitation faced by most ion trap quantum informa- 
tion experiments is scalability and to an extent, precise 
individual addressing of the ions trapped in the ground 
vibrational state of a trap. Within the LIT individual cell 
traps can be addressed with ease due to the relatively 
large separations between neighboring cells. Also, for 
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low v and low q operation the trap frequency v m is co- 
ordination number independent. In addition traps with 
the same coordination number are by definition identi- 
cal. Thus if one desires to operate with only one qubit 
per site, this system is ideally suited. 

Further the ground state of the effective harmonic os- 
cillator trap for the 3-D ion trap configuration is non- 
degenerate. However for the center of mass oscillations, 
the excited states in the symmetric operation mode of 
the trap exhibit degeneracy. Thus the usual cooling 
mechanisms (laser and sideband) essential for the pop- 
ulation of the ground vibrational state of the trap with 
the ion(s) are unchanged. In addition, specific numbers 
of ions can be loaded per site in the 3-D oscillator ground 
state. In the case of 2 ions per trap cell, proceeding with 
entangling them with available protocols using 7r/2 and 
7r pulses with appropriate detunings, or performing gate 
operations [l!^] is feasible. The one difference which arises 
due to the high symmetry is an overall rotational motion 
about the center of mass of the trapped ions which, due 
to the distance between two ions in the trap ground state 
(several /im) is very slow (ps 1Hz) compared with other 
relevant time constants in such systems. Thus entangle- 
ment of ions within a specific cell is entirely possible with 
present day technology. Applications to studies of entan- 
glement transfer, swapping [201 ] and manipulation would 
be greatly facilitated by the presence of the lattice. The 
fact that by asymmetric application of either constant 
or RF voltage the trap degeneracy can be lifted for the 
excited states offers other degrees of freedom for exploita- 
tion in the future. Thus the LIT offers a new architecture 
for quantum information processing with ions [l5j |. 

Another rich area of study with the three dimen- 
sional ion trap is the experimental few body prob- 
lem. Present day ability to load a controlled number 
of atoms/molecules/ions in the lowest energy states of 
traps allows attempts at quantitative studies with these 
very challenging systems. Linear Paul traps have long 
supported ion crystals which order in shells about the 
trap axis. The 3-D ion trap would fill in a 3 dimen- 



sionally symmetric fashion |21|, [22], [23j . Problems such 



as minimum energy configurations, phonon mode exci- 
tations and crystal relaxation, rotations, effect of trap 
symmetry and its changes can be probed. The symme- 
try of the present trap can be dramatically altered from 
a totally symmetric to completely asymmetric configu- 



ration (i.e. v„ 
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with small differences in electrode voltages. In the highly 
asymmetric configuration the system can be driven into 
chaotic behavior as evident from single, trapped ion tra- 
jectory simulations, which exhibit a dense power spec- 
trum. The ability to catch the ion signal in fluorescence 
makes trapped ions the natural choice for mapping both 
the equilibrium configurations and the dynamics of a sys- 
tem of a few, bound particles. Such studies with trapped 



ions are relevant across disciplines in physics. 

In conclusion, the concept of ion trap described here is 
versatile and suitable for many different experiments over 
and above 24| those discussed here. The trap operating 
parameters determined here are easily scaled for different 
Q/M and trap dimensions. Its primary strengths lie in 
the ease with which multiple, identical, ion traps can be 
realised and the ability to adapt the trap symmetry to 
a specific problem without additional experimental com- 
plexity. 
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